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ABSTRACT 

We spectrally fit the GeV gamma-ray flares recently-observed in the Crab Nebula by considering a 
small blob Lorentz-boosted towards us. We point out that the corresponding inverse- Compton flare 
at TeV-PeV region is more enhanced than synchrotron by a Lorentz factor square ~ T 2 , which is 
already excluding T > 200 and will be detected by future TeV-PeV observatories, CTA, Tibet AS + 
MD and LHAASO for F > 30. We also show that PeV photons emitted from the Crab Nebula are 
absorbed by Cosmic Microwave Background radiation through electron-positron pair creation. 
Subject headings: gamma rays:theory - ISM: individual: Crab Nebula 



1. INTRODUCTION 

It is well-known that the Crab Nebula is one of the 
brightest objects in the hard X-ray and gamma-ray sky. 
Because it was believed that its flux is completely steady, 
we have used it as a standard candle to calibrate detec- 
tors and instruments in those energy ranges. Since the 
Crab Nebula had already attained a position like a king 
of strong and steady sources in the high-energy gamma- 
ray sky, its impermanence must be a historic surprise. 

Quite recently AGILE (Tavani et al. 2011) and Fermi ( 
Abdo et al. 2011; Buehler 2011) reported day-timescale 
gamma-ray flares from the Crab Nebula in O(10 2 )McV- 
(D(l)GeV region, which means it is no longer stationary. 
According to the spectral fitting of the stationary compo- 
nent, the flares are most likely produced by synchrotron 
emission with an increase in the electron energy cutoff 
£max,e ~ 10 PeV and/or in magnetic field B ~ 2 mG. 
However, under the standard particle acceleration, the 
synchrotron energy loss limits the maximum synchrotron 
photons, irrespective of B, below 
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160 MeV, 
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which is violated in the flares. Possible solutions include 
the relativistic Doppler boost (e.g., Komissarov & Lyu- 
tikov 2011, Bednarek & Idee 2011, Yuan et al. 2011), the 
electric- field acceleration in the reconnection layer (e.g., 
Uzdensky et al. 2011), and the sudden concentration of 
magnetic field (e.g., Bykov et al. 2012), but there has 
been no consensus yet. 

In this Letter, we consider the relativistic model that 
a small blob is Lorentz-boosted towards us, which emits 
synchrotron radiation beyond -E™„ x (see Buehler et al. 
2011 and references therein). We stress that we can ob- 
serve the corresponding inverse-Compton flare which is 
simultaneously emitted by the same electrons existing 
in the boosted blob. Interestingly, the Lorentz factor F 
of the blob has been already constrained by the current 
TeV observations and will be further checked by the fu- 
ture TeV-PeV gamma-ray observations such as CTA, Ti- 
bet AS + MD and LHAASO, because inverse-Compton 
emission is more enhanced than synchrotron by a fac- 
tor of ~ r 2 approximately. In addition, it is remark- 



able that we must consider an absorption of PeV photons 
by Cosmic Microwave Background (CMB) radiation via 
electron-positron pair creation even for a Galactic source, 
which has not been taken into account so far. In order 
to discriminate the theoretical models and discover the 
new phenomena of the CMB absorption, the Crab Neb- 
ula is a pretty attractive experimental site for TeV-PeV 
astrophysics. 

2. STATIONARY EMISSION FROM CRAB NEBULA 

2.1. Theory and Observation 

First of all, we discuss stationary components of Crab 
Nebula emission. By assuming a broken power-law with 
an exponential cutoff for injection spectrum of primary 
electron, we parameterize it as 



dn e 
dE e = 

A e E: 
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with n e number density of electron, E e electron energy, 
B max ,e its maximum cutoff energy, s e electron spectral 
index, E c b the cooling break energy, i?ib the intrinsic 
break energy, and A e normalization. E c b is determined 
by equating the age t agc with cooling time t coo i due to 
the synchrotron energy loss. Here we assume that the 
exponent of energy on the exponential shoulder is not 
two but unity (Abdo et al. 2010). The emission below 
~ 0(1) GeV can be fitted by synchrotron radiation. The 
observational data were reported by COMPTEL (Kuiper 
et al. 2001) and Fermi (Abdo et al. 2010). We adopt val- 
ues for parameters, i age =1240 yrs, s e = 2.35, magnetic 
field B = 90/iG, the distance to the Earth d — 2.0 kpc, 
and the intrinsic breaking energy £\b = 30 GeV. Then 
the cooling energy is E c b = 1.3 TeV, and the cutoff en- 



ergy is fitted to be E n 



= 1.5 PeV. Note that the 



corresponding synchrotron cutoff energy is ~ 10 MeV, 
but the vF v peak energy is ~ 4-5 times larger than it 
because of a finite extent of the distribution. 

The emission above ~ 0(1) GeV can be fitted by 
inverse-Compton radiation due to the primary electron. 
The number density of the CMB photons is too small and 
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insufficient as target photons to fit the whoie data. Be- 
sides we also adopt the Synchrotron Self-Compton (SSC) 
process. In order to obtain target photon field for the 
SSC process, we integrate the photon number density in 
a volume where the SSC process occurs. In a one-zone 
approximation, we find 



dn 



7, target 



dr ind 2 vF u 



d 2 n. 



7, target 



V 



dtdE. 



(3) 
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with E 1 photon energy, rt 7jtarget the number density 
of photon field produced by synchrotron radiation, and 
Rssc the effective radius where the SSC process works. 
For similar parameterizations of the effective radius, 
see Atoyan & Nahapetian (1989), Atoyan & Aharonian 
(1996), and Tanaka & Takahara (2010). The observa- 
tional data in > TeV regions have been reported by 
MAGIC (Kuiper et al. 2001), HEGRA (Aharonian et al. 
2004), CELESTE (Smith et al. 2006), H.E.S.S. (Aharo- 
nian et al. 2006), VERITAS (Celik 2008, Imran et al. 
2009), CANGAROO (Tanimori et al. 1998) in addition 
to Fermi (Abdo et al. 2010), Radio and Optical observa- 
tions (Baars et al. 1997, Macfas-Perez et al. 2010). To 
simultaneously fit those data, we find the effective radius 
to be i?ssc ~ 2.0 pc. In Fig.[T]we plot the theoretical fit- 
ting and the observational data. To perform the fitting 
we use our original code which has been developed by 
one of the current authors (KK) in the series of similar 
works, e.g., Yamazaki et al. (2006). 



production is kinematically allowed with its energy ex- 
ceeding threshold E lb > m 2 c A /(7Ej). The attenuation 
length Q is given by 



fatten 
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and the pair-production cross section through 7 + 7 
e + + e~ is given by 



&pair(E^ : E^ h , /l) 
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with /3 = y/l - Am 2 e c 4 /s and s = 2£ 7 £ 7b (l - /t). When 
we consider a 10 3 TeV (1 PeV) photon, the threshold 
energy of the target photon for the pair production be- 
comes the order of O(10~ 3 ) eV at which the CMB photon 
dominates along the line of sight between Crab Nebula 
and the solar system. In Fig. [5] we plot the attenuation 
length as a function of the photon energy in eV. Remark- 
ably the attenuation length can be down to 4 - 5 kpc at 
E 1 ~ (2 - 3) PeV. 
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Fig. 1. — Spectrum fitted to the stationary component of radia- 
tion from Crab Nebula. The solid line shows the total spectrum. 
The dashed and the dotted lines represent the synchrotron and the 
inverse-Compton emissions, respectively. The upper and the lower 
curves of inverse-Compton process at TeV regions are for scattering 
off the synchrotron and CMB photons, respectively. Observational 
data are plotted as points with their error bars. We also show 
sensitivities of future projects such as CTA (Actis et al. 2011) 
LHAASO (Cao et al. 2010) and Tibet AS + MD (Takita 2011) for 
50 hours measurements denoted by (50h). 

2.2. PeV gamma-ray absorption by CMB 

Photon is absorbed if there is a sufficient number 
of background photons and the electron-positron pair 




Fig. 2. — Gamma-ray horizon as a function of energy for inci- 
dent photon. In this energy region, the electron-positron pair pro- 
duction through the scattering off the background CMB photon 
dominates the energy-loss rate. 

Because the CMB radiation is ideally homogeneous 
and isotropic, we can simply calculate the absorption fac- 
tor by using exp[— d/ L at tcn{E.y)} with d — 2.0 kpc (Abdo 
et al. 2010). Multiplying by this absorption factor, we 
obtain the observed spectrum. In Figj3] we plot spectra 
with and without this type of absorption by the CMB 
radiation. In addition, we also plot sensitivities of fu- 
ture projects such as LHAASO (50 h and 100 h), CTA 
(50 h and 500 h) and Tibet AS + MD (50 h and 500 
h) with their observation time written in round brack- 
ets, respectively. In Fig. [4] we compare the observational 

1 The attenuation length is equal to the interaction length for 
our parameters and the Calactic magnetic field. 
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Fig. 3. — Gamma-ray spectra with and without the CMB ab- 
sorption. We also plotted sensitivities of future projects such as 
LHAASO (50 h and 100 h), CTA (50 h and 500 h) and Tibet AS 
+ MD (50 h and 500 h) with their observation time, respectively. 
Note that the vertical axis means vF v times E^. 

sensitivities with the net absorbed component, which is 
the difference between spectra with and without the ab- 
sorption. From this figure, we find that the absorption 
by CMB radiation will be observed by LHAASO with 
its observation time of 100 hours. In turn, this means 
that we have to consider this new type of absorption by 
the CMB radiation whenever we observe the PeV photon 
from Crab Nebula. As far as we know, we point out this 
detectability of the PeV photon absorption by the future 
telescopes for the first time. 



3. FITTING TO FLARE COMPONENT 

Recently some observations reported that Crab Neb- 
ula is no longer stationary with flares at around 0(1) 
GeV (Tavani et al. 2011, Buehler 2011, Buehler et al. 
2011). Fig. [5] shows these data points. The duration 
A£ D bs is typically the order of one day. In this section 
we discuss how we can explain these flares in terms of 
synchrotron emission by accelerated electrons. We con- 
sider Lorentz-boosted blob models in which a small blob 
is boosted with a Lorentz factor T and an off-axis view- 
ing angle 9. In addition, as will be discussed later, it 
should be natural to assume accelerated electrons and 
magnetic field in the blob. Electrons emit synchrotron 
radiation by using the local magnetic field in the rest 
frame of the blob. This model can produce 0(1) GeV 
synchrotron photon in the observer frame unlike the non- 
relativistic models where the synchrotron photon energy 
cannot exceed E™ x ~ 10 2 MeV in equation (Q}. -E™„ x is 
independent of B because -E ma x,e is limited by balancing 
synchrotron cooling time with acceleration time. Here 
we do not specify the origin of the blob, which could be 
the pulsar wind or the shock at the knot of the inner 
nebula. 

Importantly, inverse- Compton photon is also emit- 
ted simultaneously, which is produced by scattering off 
both the boosted CMB and synchrotron radiation com- 
ing from outside the blob. Both of synchrotron and 
inverse- Compton radiation are boosted-back to the ob- 
server frame, which give larger energy and sizably-larger 
fluxes even if the size of the blob is tiny. 

Concretely, the photon energy flux vF v (in unit of erg 
cm -2 s _1 sr _1 ) emitted in the rest frame of the blob can 
be boosted in the observer frame by a following scaling, 



vF^E-,) -> S l 'uF v {R 1 /S), 
where the Doppler factor is represented by 



(7) 




Fig. 4. — Same as Fig.[3]but for the difference between two lines 
(with and without absorption). Note that the vertical axis means 
Ej times A(uF^) in linear scale. 



r(l-j9cos0)' 



(8) 



which behaves like ~ T for 9 < 1/r, and ~ 1/r in the 
other limit for off-axis cases, 9 ^> 1/r. Indeed, this shift 
of energy for synchrotron emission can fit the change of 
cutoff scale for the flare component. In addition, for 
inverse-Compton processes in the rest frame of the blob, 
the target photon is boosted from the observer frame, 
and its distribution is modified to be 



dE. 



dn~, 



{E 1 



TdE 1 ^{E 1 /T). 



^ dE^ ^ ^ dEry 



(9) 



Therefore the inverse-Compton power is more enhanced 
than the synchrotron power by a factor of T 2 for the 
Thomson limit, 



iyF v \ 



IC 



IC 



iyEy )syn 



(10) 



Considering the Klein-Nishina effect, the enhancement 
could be smaller than T 2 . 

Next we discuss our choice of model parameters. If the 
maximum energy of electron in the blob rest frame is lim- 
ited by the synchrotron cooling, the bulk Lorentz factor 
should be < 10 to explain the energy shift of the maxi- 
mum energy -E™f x shown in equation ^ at the flare. In 
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that case, the synchrotron cooling time should be shorter 
than the variability timescale, which gives 



EUe > 170 TeV ( 



_1 fAt 



bs 



V 8h 



where E' m&x e is the maximum energy of electrons, and 
B 1 is magnetic field in the rest frame of the blob, respec- 
tively. 

The Lorentz factor can be larger than > 10 because the 
maximum energy of electrons is not necessarily limited 
only by cooling, but by the blob size (e.g., Ohira et al. 
2011). In this latter case, the Larmor radius of electron 
with the maximum energy would be comparable to the 
blob size, ~ 5cAt t, s in the blob rest frame, which gives 



EL 



= 790 TeV 



D' 



At bs 



(12) 



Equating (jTTJ) with (|12p , we obtain a threshold magnetic 
field in the blob rest frame, 



(13) 



For B' > B' th , the maximum energy of electron is lim- 
ited by the synchrotron cooling. However, B' th is much 
larger than that expected from the standard model of 
the steady-state Crab Nebula (e.g. Kennel & Coroniti 
1984). Therefore, we also consider the case, B' < B' th , 
where the maximum energy of electron is limited by the 
blob size. Because the observed energy of synchrotron 
photons for a monoenergetic electron is written as 



£ syn = 95 MeV ( |() 



EL, 



500 TeV 



B' 
3 mG 



, (14) 



by removing E' max e from Equation (TT2"j) . the magnetic 
field in the blob rest frame is obtained as 



B' = 2.2 mG 



E. 



syn 



10 2 MeV 



1/3 



At 



obs 
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Then the condition B' < B' th gives lower bound on the 
Doppler factor 



5 > 22 



E S yn 



10 2 MeV 



(16) 



By removing 6 in (|12l) with (|14p . we get the maximum 
energy of electron in the blob rest frame, 
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10 2 MeV 
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(17) 



which depends only on observables. In this size-limit 
case, the variability timescale should be determined by 
the dynamics of the blob, since the cooling timescale is 
longer than that. This may be favorable for explaining 
the comparable timescales of rise and decay in the ob- 
served flares. 

In Fig [5] we plot theoretical calculations of the spec- 
trum fitted to the flare component, adopting a scaling 
law for magnetic field B' = 220/iG(£/10 2 )~ 1 in equation 
(|15p and the maximum energy of electron obtained in 



equation (fTT)) . The thick solid line shows a result with 
r = 5 = 10 2 and B = 223 fxG. We call this set of pa- 
rameters "fiducial model". Because the electrons in the 
blob cannot be cooled in such a short time, the cooling 
break does not appear in GeV region. Then the spec- 
tral index of electron could be the same as the initial 
value, n s ~2.35. Consequently the photon index of syn- 
chrotron radiation (F v jv) at around O(10 2 ) MeV be- 
comes (n s + 1)/2 ~ 1.6, which is consistent with the ob- 
servation (1.27±0.12) reported by Buehler (2011). Below 
O(10) MeV, the synchrotron flare component is smaller 
than the stationary one. We also plot the high and low 
r models with r = 5 = 700 and r = 5 = 30 by using the 
same scaling law for the magnetic field. 

On the other hand, in TeV-PeV region, it is remark- 
able that the inverse- Compton radiation of the flare com- 
ponent can exceed the stationary one significantly at 
E 7 > O(10) TeV for the fiducial model. That is also due 
in part to that the inverse-Compton power is enhanced 
more than the synchrotron power because of the boosted 
target photon distribution by a factor of T 2 shown in 
PH|) . 0. From Figure [5J we find that higher Lorentz fac- 
tor than r > 700 has been already excluded by the cur- 
rent GeV observations and T > 200 by TeV observations 
(Mariotti et al. 2010; Ong et al. 2010). 

The future gamma-ray observatories such as CTA, Ti- 
bet AS + MD and LHAASO will be able to probe T down 
to ~ 30 by 8 hours observations (by one night). The sen- 
sitivities are shown in Fig. [5] by conservatively linearly 
scaling them from 50 hours to 8 hours. Even if the flare 
flux is less than the stationary one, we can detect it down 
to the statistical error of photons. Note that the ratio of 
inverse Compton to synchrotron power depends only on 
r in equation (|10p . not on 8, so that we can constraint 
r even for an off-axis event. Note also that the cooling 
effects can be seen a little in > PeV region for T < 30. 

To be more conservative, we also show similar cal- 
culations with adding an artificial low-energy cutoff for 
the electron spectrum where we took zero flux for E e < 
E e ,cut = 10 2 TeV in the blob rest frame, which is repre- 
sented by the shallower thick solid lines [this may explain 
the X-ray feature (Tavani et al. 2011)]. Even in this case, 
we can probe T down to ~ 700. 

Therefore we can discern the model from the nonrela- 
tivistic ones for the flares at 0(1) GeV by observing the 
inverse-Compton radiation in the O(10)TeV-O(l)PeV 
energies. 

Here it should be meaningful to check an energy ra- 
tio of electron to magnetic field in the blob. The flare 
luminosity in the blob rest frame, L', is given by 



II = ATtd 2 {vF v ) ohs 8- 



(18) 



Assuming that this luminosity is produced by syn- 
chrotron radiation of electrons with a typical maximum 
energy, the luminosity is written as 



L' = K(E' m 



, 4 / B' 2 



-EL 



m e c* 



(19) 



2 When we use the scaling B' <x <5 _1 in equation H15I I in order 
to fit the synchrotron component in Ge V re gion, the magnetic en- 
ergy density is reduced by 8~ 2 . With mot , the inverse-Compton 
component is enhanced by a factor of T 2 S 2 in total. 
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where N' c {E' max e ) is the number of electrons with the 
maximum energy in the blob rest frame. From Equa- 
tions ([T5)l. (fT7| . ([18]) and ([TU), we obtain the total energy 



of electrons with E' av . 



u ' e = N' e (E' m )E' m 



1.3 x 10 37 erg 



10 2 MeV 



-7/6 



8 h 



(20) 
1) 



On the other hand, from equation (|15|) . the total energy 
of the magnetic field in the blob rest frame is given by 



B' 2 4 3 
— — X -7T (<ScAt obs ) 

OTT O 



= 4.8 x 10 4 
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(23) 



for the size-limit case B' < B' th - Then, we find the ratio 
of the electron energy to the magnetic energy to be 



EL 
u' B 



= 2.6x10" 



10 2 MeV 



-11/6 



At 



obs 



-2/3 



8 h 

This ratio can increase by ~ (E maXi< ,/ Eib) ns 2 for 2 < 
n s < 3 if we include the low energy electrons. Even 
in this case, the boosted blob might be magnetically 
dominant. Since the blob position, TScAt h s ~ 8.6 x 
10 l4 T<5 cm, may be also much smaller than the radius of 
the termination shock, ~3x 10 17 cm, the blob might be 
produced in the pulsar wind before the magnetic energy 
is converted to the bulk kinetic energy. Alternatively, the 
observed flares might be off-axis ones with small 5(~ 1). 
In this case we predict larger flares than ever detected. 
Or a flare may consist of many pulses with At b s < 8 h, 
e.g., U' e ~ U' B for At bs ~ 33 msec (pulsar period). It 
may be the reason for the scarcity of flares that a flare 
needs many pulses. 

If we require that the blob energy TU' B is 
less than the spin-down energy during the flare 
EAt b s 5/r in the observer frame, we have T 2 < 
300(£ syn /10 2 MeV)- 2 / 3 (At obs /8h)- 2 / 3 for the size-limit 
case B < B' th . Thus, if we find T > 30, we also imply 
that a flare consists of many pulses with Ai b s < 8 h. 

As was mentioned in Introduction, so far there has 
been no consensus of the theoretical models. In the 
model with only increasing the maximum energy of elec- 
tron such as by the electric acceleration, there is surely 
an excess in PeV region for inverse-Compton flare com- 
ponent. In order to detect this excess by LHAASO, how- 
ever, we approximately need a few tens of hours for the 
observation time, which is longer than the typical dura- 
tion of the flare. In the model with only increasing the 
magnetic field such as by rapid compression, the inverse- 
Compton flare is highly suppressed. 

4. SUMMARY AND CONCLUSION 

In order to explain the origin of the GeV flare in Crab 
Nebula, we have studied models in which a small blob 
is boosted, e.g., with a Lorentz factor T > 30, and 
emits synchrotron photon higher than the maximum syn- 
chrotron energy E™^ shown in equation (TTJ). We have 
also discussed possibilities that we will discern the model 
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Fig. 5. — Spectrum fitted to the flare component of radiation 
from Crab Nebula. We adopted observational data for the flare 
reported by Fermi (Buchlcr 2011, Buehler et al. 2011) which is 
denoted by Fermi2011. The thick solid line shows the theoretical 
prediction with r = 6 = 10 2 , The thin solid lines are the cases for 
r = <5 = 30 and F = 5 = 700, respectively. The corresponding 
models with low-energy cutoff (B e ,cut = 10 2 TeV) are represented 
by the shallower solid lines. The meanings of the other lines and 
the observational data are the same as those in Fig. [l] 

from the others such as nonrelativistic models, by ob- 
serving the corresponding inverse-Compton flare compo- 
nent. We have pointed out that the inverse-Compton 
flare can appear in > 0(1) TeV region accompanied with 
the GeV flare in this kind of the boosted blob models with 
large Lorentz factor because the inverse-Compton power 
is more boosted than the synchrotron power by ~ T 2 . 
High r models have been already excluded for T > 200 
by the current TeV observations and will be further down 
to r ~ 30 by the future TeV-PeV observatories, such 
as CTA, Tibet AS + MD or LHAASO. In addition, by 
considering this enhancement in the TeV-PeV region, in 
near future we may observe "orphan TeV-flares" , which 
do not have even a GeV flare. 

Even for the stationary component of Crab Nebula, 
we have also pointed out for the first time that the ab- 
sorption of PeV photons by CMB radiation through pair 
creation 7 + 7cmb — > e + + e~ is important. We must 
consider this effect whenever we fit the spectrum of Crab 
Nebula in the C(l)PeV regions. 

It is notable that we will be able to accomplish those 
studies for observation of Crab Nebula at 0(l)TeV- 
0(l)PeV energies by using the future gamma-ray tele- 
scopes such as CTA, Tibet AS + MD or LHAASO. We 
hope the earliest possible completions of this kind of new 
gamma-ray telescopes. 
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